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ABSTRACT 

New radial velocity measurements for previously known and newly confirmed globular clusters in 
the nearby massive galaxy NGC 5128 are presented. We have obtained spectroscopy from LDSS- 
2/Magellan, VIMOS/VLT, and Hydra/CTIO from which we have measured the radial velocities of 
218 known, and identified 155 new, globular clusters. The current sample of confirmed globular clusters 
in NGC 5128 is now 605 with 564 of these having radial velocity measurements, the second largest 
kinematic database for any galaxy. We have performed a new kinematic analysis of the globular 
cluster system that extends out to 45' in galactocentric radius. We have examined the systemic 
velocity, projected rotation amplitude and axis, and the projected velocity dispersion of the globular 
clusters as functions of galactocentric distance and metallicity. Our results indicate that the metal- 
poor globular clusters have a very mild rotation signature of 26 ± 15 km s . The metal-rich globular 
clusters are rotating with a higher, though still small signature of 43 ± 15 km s _1 around the isophotal 
major axis of NGC 5128 within 15'. Their velocity dispersions are consistent within the uncertainties 
and the profiles appear flat or declining within 20'. We note the small sample of metal-rich globular 
clusters with ages less than 5 Gyr in the literature appear to have different kinematic properties than 
the old, metal-rich globular cluster subpopulation. The mass and mass-to-light ratios have also been 
estimated using the globular clusters as tracer particles for NGC 5128. Out to a distance of 20', we 
have obtained a mass of (5.9 ±2.0) x 10 11 M Q and a mass-to-light ratio in the B-band of 16 M q /Lb0- 
Combined with previous work on the ages and metallicities of its globular clusters, as well as properties 
of its stellar halo, our findings suggest NGC 5128 formed via hierarchical merging over other methods 
of formation, such as major merging at late times. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 5128) — galaxies: 
kinematics and dynamics — galaxies: star clusters — globular clusters: general 



1. INTRODUCTION 



The modern idea of galaxy formation involves 

hierarchic ally merging smaller galax ies or proto- 

galaxie s dToomrd Il977t iWhite fc Reel [19781: iPeeblesI 



1980; IWhite fc Frenkl I 19911: iKauffmann et all 11995.. 
Baughetall 119981: ICole et al.l 120001: ISomerville et al I 
20011 ) that contain both baryonic and dark matter. The 
merging of these smaller clumps into larger systems re- 



matiasKomcz@unab.cl 
harrisSSphysics.mcmaster.ca 
dgeisler@astro-udec.cl 



glharrisSS astro, uwaterloo.ca 

Current Address: Department of Physics & Astronomy, Uni- 
versity of British Columbia, Vancouver BC V6T 1Z1, Canada; 
kwoodley@phas.ubc.ca 



suits in a bound system embedded in a large dark matter 
halo. 

Globular clusters (GCs) are great tracers of the 
formation and evolution of their host galaxy as there is 
significa nt evidence that they trace episodes of star for- 
mation (fHoltzman et al.lll9~9"2t ISchweizer fc Seitzerl 19931 : 



Whit more et all 119931: 



lweizer et al 



Zenf et al l 119951: ISc 
19971 iCarlson et al.l"1l998t 
Whitmore et all 119991: IZepf et al 



Whitmore fc Schweizerl 
~li996l: 



199E 



Miller et all 
Schweizer fc Seitzerl 11 998: 
19991: iChien et ail 



20071: IGoudfrooii et all 120071: iTrancho et al.l 120071 

among others), and generally have ages > 10 Gyr, 
making t hese some of the oldest structures in the 
Universe (Kisslcr- Patig et all 119981 : iForbes et all 120011 : 
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Woodley et al. 



Puzia et alJ l 200llHempel et al.ll200l Kundu et al.ll2005t 



Puzia et al.l 120051: IStrader et all 120051: iSharina et al I 



20061: IBeaslev et all 120081: iPuzia fc Sharinal 120081: 
Woodlev et al.1 12010L as examples). GCs are also an 



excellent way to study massive galaxies because they 
are luminous, compact and can be found in large 
quantities out to large galactocentric radii in relatively 
spherical distributions. They are also expected to 
retain a memory of their formation conditions (e.g. 
IPuzia et al.ll2006f ). so by understanding their formation 
we may be able to constrain formation scenarios of 
their host galaxy. Finally, each cluster has a unique age 
and abundance which are relatively straightforward to 
determine, providing us with the key data we need to 
trace their formation and chemical evolution. 

The color bimodality, consisting of GCs having 
two distinct peaks in their broadband color his- 
tograms, corresponding to red and blue GCs, is a 
common feat ure in all e arly-t y pe galaxies such a s 
NGC 4472 (ICohen et all [20031: iStrader et all l200l. 



JL 

M87 ([Kundu fc Zepfll2007l). a nd NGC 5128 (IFeng et al.l 

~7u 



l2004bT IBeaslev et al.1 2008 IWoodlev et al .1 120 ion . as 
well as late-type galaxies, such as the Milky Way and 
M31 (IZinnl 119851: lElson fc Walterbosl Il988t lArmandroffl 
1981 iBrodie fc Huchral 119901: iZinnl 119901: iHarrisI 119911: 
Ashman fc Zepj I1998D . The separate populations of 



GCs have led to the idea that they form via different 
mechanisms /condition s and perhaps at different times. 
lAshman fc Zepj ([19921 ) proposed that red, or metal-rich, 
GCs can result from the merging events of gaseous disk 
galaxies. This type of formation suggests an age spread 
in the metal-ri ch GCs as disk merg ers ca n continue to 
low re dshifts. IBeaslev et al.l ([2002) and IStrader et al.l 
(2005) have proposed an interesting formation scenario 
that com bines two scenarios proposed early: multiphase 
collapse ([Forbes et al.lll997T ) which suggested that most 
(if not all) GCs form in their parent galaxy, but that 
the galaxy undergoes more than one phase of collapse 
and hence star format i on: a nd the accretion scenario 
([Cote et all I1998L I2000L 120021 ) suggesting that blue, or 
metal-poor, GCs are accreted into the massive galaxy 
via merging of surrounding satellite galaxies. Within this 
scenario, the metal-poor GCs could form in small proto- 
ga lactic clouds i n the early universe. As first suggested 
by [Santos (200f), their formation could have been trun- 
cated by the reionization process. It is also possible that 
the star formation was truncated by feedback processes 
such as the removal of gas from supernovae drive n winds 
(e.g. IDekel &TWodi\2Qm Ikauffmann et al.ll2003D . These 
protogalaxies would merge and form a massive early-type 
galaxy with the field stars and metal-rich GCs forming 
in the second collapse of star formation. 

In any of the above, it is clear that age is a dis- 
tinguishing factor. However, the metallicity bimodal- 
ity also suggests the blue and red components may 
have different kinematical signatures. Kinematics can 
be used to test whether red and blue GCs share sig- 
natures of rotation and velocity dispersion, and how 
those signatures can vary with galactocentric distance. 
We already see in many elliptical galaxies that the red 
GCs are mor e centr a lly concentra t ed than the blue (e.g 



question whether their kinematics also change with ra- 
dius. 

NGC 5128 (Cent aurus A) is a mod erately luminous 
(M B = -21.1 mag. lDufour et al|[l979T) elliptical galaxy 
in the nearby Centaurus group, at a distance of 3.8 ±0.1 
Mpc ([Harris et al.l 120091) . It has evidence for a rela- 
tively recen t merger ev ent such as a warped disk and 
faint shells ([Ma lin 1978 ), r e cent s tar formation (jGrahaml 
[19981: IReikuba et all 120011. [200l . and a small blue el- 
liptical arc that appears to be a tidally disrupted stel- 
lar stream, which could be a recently accreted satel- 
lite galaxy ([Peng et al.l 120021) . The close proximity of 
this galaxy allows a detailed study of its stellar popu- 
lation, including the GCs, that is not yet possible for 
other giant elliptical galaxies. However, this close prox- 
imity also means the galaxy is spread across a very 
wide area in the projected sky. NGC 5128 also has 
a low galactic latitude of b = 19° so its GCs need to 
be securely confirmed via their radial velocities, which 
are generally distinct from the two forms of contam- 
ination, Milky Way foreground stars and background 
galaxies. Currently, over 400 GCs have been confirmed 
with this methodolog y in t h e past (Ivan den Bergh et all 
19811: iHesser et all 119841 119861 : Harris et al l Il992t 



Peng et al.l l2004bl: IWoodlev et all [20051: IReikuba et al 



Geisler et al.l [l996t IForbes et al.l [l998t iRhode fc Ze 
20011: iDirsch et all 120031: iForbes et all 120041: iPeng et al 



2004b; iTamura et all 120061: IWoodlev et alll2007D . so we 



20071; IBeaslev et all 120081: IWoodlev et all 120101 ) . with 
a large number of these compiled in the catalog of 
IWoodlev et all j2003). 

In this paper, we confirm a number of new GCs as 
well as remeasure previously known GCs in NGC 5128. 
The observations and data reduction are presented in 
Section followed by our presentation of radial velocity 
measurements in Section [3] In Section [4l we present 
a new kinematic analysis of the GCs in NGC 5128 as 
well as estimate a new mass and mass-to-light ratio of 
NGC 5128. We then discuss our findings in Section [5] 
and conclude in Section [BJ 

2. OBSERVATIONS AND DATA REDUCTION 

Our spectroscopy involves three independent observ- 
ing runs, each with a primary objective to obtain radial 
velocity measurements of GC candidates to confirm new 
membership. Table [I] lists the center of field location in 
R. A. and Dec! (in J2000 coordinates) and the total 
exposure times for each field. 

2.1. LDSS-2/Baade Dataset 

We obtained multi-object slit spectroscopy with 
the now decommissioned Low Dispersion Survey 
Spectrograph-2 (LDSS-2) on the Baade/Magellan 6.5 m 
telescope at the Las Campanas Observatory in Chile. We 
observed 12 fields scattered around the central regions of 
NGC 5128 on May 9-10, 2002 (PI: D. Geisler) with the 
6' field of view of LDSS-2 (at the distance of NGC 5128, 
1' is approximately 1.1 kpc). Our slit sizes easily encom- 
passed both our target objects as well as the sky imme- 
diately surrounding each object. A gain of 1.0 e~/ADU 
and a readout noise of 5.0 e~ were used in this study. 

Our general reduction was performed with the Image 
Reduction and Analysis Facility (IRAF0 software. Each 

2 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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of our target fields was taken with exposure times be- 
tween 1000-1800 seconds. We obtained and combined 
zero (bias) frames (zerocombine) for each night and re- 
moved the bias from our science exposures. The over- 
scan region of the CCD was also removed from our sci- 
ence target exposures (ccdproc), followed by the removal 
of cosmic rays with the European Southern Observatory- 
Munich Image Data Analysis System (ESO-MIDAS) task 
filter/cosmic. The flat fields were combined for each field 
(flatcombine in IRAF) and applied to our science fields 
(ccdproc). Within each slit we fit, traced, and extracted 
each object and sky spectrum (apall). Comparison arcs, 
taken adjacent to field exposures, were used to gener- 
ate a transformation from pixel to wavelength for each 
of our target objects (identify) and applied with dispcor. 
We used the neutral oxygen 5 5 77 A sky line to correct for 
small shifts (typically less than 0.5A) in our spectroscopy 
(specshift) and then subtracted the sky from the targets 
(sarith) and then cleaned the spectra (lineclean) using a 
low rejection of 8 and a high rejection of 2.5 below and 
above the residual sigma with a spline3 function. Each 
field was observed two times and then sum combined 
(scombine). Our radial velocities were then measured 
with fxcor using GC0106 as a template spectrum. This 
object is a bright cluster with previously well determined 
radial velocity from the literature. 

2.2. VIMOS MOS/VLT Dataset 

We obtained spectroscopy with the VIMOS MOS/VLT 
8 m telescope in 2007 (PI: M. Gomez, Program ID 079.D- 
0539), on Paranal, Chile. We obtained 5 fields located 
around the inner regions of NGC 5128, each with 4 CCD 
quadrants of 7' x 8' separated by 2' gaps. We targetted 
166 previously known GCs and the remaining slits were 
placed on candidate GCs. Our slit width was 1.0", em- 
compassing both our target objects as well as surround- 
ing sky. 

Our target data fields were taken in ~ 1000 sec- 
ond exposures with the high-resolution (HR) blue grism 
with no filter and approximate wavelength coverage of 
4100 - 6300A. The spectral resolution was 2050 and the 
dispersion 0.51 A/px for a 1" slit. We observed with 
a binning of 1 by 1 and a gain of 1.96 e~/ADU and a 
readout noise of 3.1 e~. 

Our general reduction performed with IRAF involved 
combining (zerocombine) and removing our zero (bias) 
exposures from the science targets along with subtract- 
ing the overscan region (ccdproc). We then subtracted 
the cosmic rays in our target fields with the ESO-MIDAS 
task filter/cosmic. Following the cosmic ray removal, we 
returned to IRAF to identify and fit the sky regions as 
well as trace and fit each target spectrum with a legen- 
dre function of order 5. We then extracted and removed 
the sky spectrum from our target spectrum (apalt). A 
comparison arc, taken with each target field, was used to 
identify lines (identify) from which we applied the dis- 
persion correction to the data (dispcor). The different 
exposures for each field were combined (scombine) and 
the radial velocitie s were measured (f xcor) using the tem- 
plate GC 158-213 (jPuzia et al.ll2002h . 

2.3. Hydra/CTIO Dataset 



We obtained spectroscopy from the Hydra Spectro- 
graph on the Cerro Tololo Interamerican Observatory 
(CTIO) Blanco 4 m telescope, in Chile (PI: W. Harris, 
Program ID 2007A-0070). We observed from April 7 to 
April 10, 2007 and had 2.5 usable nights of data out of 
4, with lost time due almost entirely to poor weather 
conditions. We observed 4 target fields, each with a 40 
' field. The fiber spectrograph has 2 " diameter fibers, 
easily capturing the light of a GC at the distance of NGC 
5128. We placed available fibers on both candidate and 
known GCs as well as approximately 25 fibers placed on 
random sky positions in each field. Our sky fibers were 
scattered throughout the full radial range and within the 
full concentric range of the field. 

All of our data was obtained in the red regime with a 
range of ~ 2300A centered on 7600A, capturing Ha and 
the two bluest of the Calcium triplet lines. We used the 
KPGLF grating with a dispersion of 0.57A/px and no 
filter. Our binning was 1 x2 (2048 x 4096 pixels) and we 
used a gain setting of 0.84 e~/ADU and a readout noise 
of 3 e~ . 

We reduced our data with the IRAF software. We re- 
moved the combined zerocombine bias in IRAF as well as 
subtracting the overscan region of the CCD from the ex- 
posures (ccdproc). The cosmic rays were r emoved using 
the L.A. Cosmic task (Ivan Dokkum 200l]). We created 
our flat field using Milky flats taken on the first day of 
observations by combining our 5 exposures (each of 300 
seconds) using flatcombine. The task dohydra was then 
run on each target exposure to define the apertures, fit 
a flat field to our combined milky flat, as well as define a 
pixel-to- wavelength transformation with our comparison 
arc spectrum, which was then applied to our target data. 
Next, we combined our twilight exposures for each night 
(flatcombine) and corrected for the fibre throughput dif- 
ferences. During our pixel-to-wavelength correction, we 
also trimmed the exposures to include 6500 - 8740A with 
2000 pixels. We then extracted the spectra from multi- 
to one-dimensional fits file format with scopy. 

In order to remove the sky contribution, we scaled the 
averaged sky spectrum for all fibers in each field to match 
the intensity of the spectrum of each individual object, 
which was then subtracted from the target spectrum us- 
ing sarith. We removed the continuum from each expo- 
sure (continuum) with a 10th order legendre function, 
rejecting the 6 lowest and 3 highest pixels. Each object 
was then combined with the same target from all expo- 
sures with scombine by taking the average. 

We measured radial velocities of our target objects by 
creating a template GC using two bright GCs with previ- 
ously measured radial velocity measurements from high- 
resolution spectroscopy. We chose GC0330 and GC0365 
from field 1 with weighted velocity measurements of 
673±1 km s _1 and 594±1 km s -1 , respectively. We com- 
bined the two fully reduced, sky subtracted spectra by 
averaging with scombine after each spectrum had been 
doppler corrected to km s _1 . With our template in 
hand, we used fxcor to cross-correlate each object with 
our template M31 cluster GC 158-213. 

3. GLOBULAR CLUSTER SYSTEM OF NGC 5128 
3.1. Newly Confirmed Globular Clusters 



4 



Woodley et al. 



We have remeasured radial velocities for 218 pre- 
viously known GCs, listed in Table [2j This table 
lists the known GC ID, the R.A. and Decl. (J2000), 
and the velocity measurements presented here. It 
also lists the new weighted veloc ity measurement of 
all GCs in the past literature (Ivan den Bergh et al.1 



l98l IHesser et all fl98i 



(jva 

1986; JHarris et al 
"all " 



1992; 



2004d: IWoodlev et al.l 120051: iReikuba et al.1 



Pone et al. . _ 

2007tlBeaslev et al.ll2008t IWoodlev et al.ll2010D including 



the new measurements of this study. We have remea- 
sured 64 GCs with LDSS-2, 154 GCs with VIMOS, and 
57 GCs with Hydra. Least square fits are shown for the 
comparison of our measured velocities with known val- 
ues in the literature in Figure [TJ The comparisons for 
nearly all GCs are well within the measurement uncer- 
tainties. In the rare cases of outlying measurements, the 
most obvious being GC0445 in the Hydra study, there 
is only one previous radial velocity measurement from 
IWoodlev et al.l (j2010fl . but both measurements are still 
within the bounds of accepted radial velocity confirma- 
tion of 200 — 1000 km s _1 . We have higher confidence 
in the radial velocity measured with CMOS because the 
spectrum had a higher S/N than the Hydra data. 

Table[3]lists our 155 newly confirmed GCs from our ra- 
dial velocity measurements. Table[3]lists the GC ID, the 
R.A. and Decl. (J2000) the C, M, and Ti magnitudes 
from lHarris et al.l ((2004al ) , and the measured radial ve- 
locity. We identified 29 new GCs with LDSS-2 (GC0416- 
GC0444), 80 new GCfOwith VIMOS (GC0480-GC0559), 
and 46 new GCs with Hydra (GC0560-GC0605). The fi- 
nal column in the table, lists the weighted average of all 
radial velocity measurements of these new GCs. 

We compare all multiple measurements of the same GC 
within each radial velocity study to test for consistent 
measurements. A comparison fit is shown in Figure [5] for 
the 4, 1, and 17 multiple measurements for the LDSS- 
2, VIMOS, and Hydra studies, respectively. An average 
uncertainty for the studies is also indicated in the figure. 
We find the multiple measurements match a 1:1 fit within 
the measurement uncertainties of each study. The mean 
rms difference between multiple measurements are 64 km 
s- 1 for LDSS-2, 22 km s" 1 for Hydra, and 35 km s" 1 
for VIMOS (recall there is only 1 GC measured in the 
VIMOS dataset that was measured in multiple fields). 

Of our measured GC candidates there are a small num- 
ber with velocities between 150 — 250 km s _1 which may 
be either GCs or foreground Milky Way stars. This ve- 
locity domain is the major region of contamination of the 
GC catalog in NGC 5128. To distinguish between these 
two objects, we measured their structural parameters us- 
ing our superb 1.2 x 1.2 deg 2 Magellan/Inamori Mag- 
ellan Areal Camera and Spectrograph (IMACS) images 
(R filter) taken in 0.5" seeing. The structural param- 
eters were measured with ISHAPE (|Larsen "et~aLllT99l 
120011 ) which convolves the stellar point s pread (PSF) 
function with an analytical King profile (|Kingl 119621 ) 
and compares the result with the input candidate im- 
age achieving a best match. The structural parameters, 
measured from the models, that we used are the core 

3 Four of these newly identified GCs may be stars based on their 
structural paramter measurements, but have velocities within the 
accepted range. These objects, GC0485, GC0508, GC0523, and 
GC0526, are included here as new GCs but with caution. 



radius, r c , the tidal radius, r t , the concentration pa- 
rameter, c = r t /r c , and ellipticity. The half-light radii 
can also be obtained from the transformation, r e /r c ~ 
0.547 c - 486 which is good to ±2% for c > 4 (|Larsen et al.l 
2001). which is satisfied for GCs in NGC 5128 (see 



Gomez fc W oodley 2007]). We compared the measured 



structural parame ters of the GC candidates to normal 
GCs in NGC 5128 Ipijris et al.ll2006UG6mez et alJ l2006l 
iGomez fc Woodlevll2007t iMclaughin et al.ll2008| ). as well 
as to a delta distribution light profile, which describes 
the PSF of a star. Based on their structural parame- 
ters, we classified a small handful of candidates as GCs 
with low velocity measurements. With our preliminary 
look at the structural parameters of GCs in NGC 5128, 
we note that five GCs are well fit by a delta distribu- 
tion light profile. Four of these are newly identified here 
(GC0485, GC0508, GC0523, and GC0526), while the 
other is GC0472. A closer look at the structural pa- 
rameters of these objects is underway, however, we have 
included them here because they satisfy the velocity cri- 
teria for the GCs in NGC 5128. We have determined the 
kinematics below including and excluding these 5 GCs 
and found that there is no kinematic difference between 
the two solution sets within the uncertainties. The re- 
sults presented thus include the objects. 

3.2. Radial Velocity Distributions 

We have divided our globular cluster system (GCS) 
into either metal-poor or metal-rich GCs based on 
their color information. We calculated a metallicity 
[Fe/H]c_Ti obtained fro m a transformati o n of dered- 
dened (C-Ti) derived by lHarris fc Harris! ((2002) cali- 
brated through the Milky Way GC data. The trans- 
formation, calibrated in the range of (C-Ti) = 0.9 to 
1.9 mag and [Fe/H]= -2.3 to -0.1, is 

[Fe/H] C - Tl = -6.037x(l-0.82x(C-T 1 ) o ) + (0.162x(C-T 1 ) 2 ) 

(1) 

where 

(C-Ti) = G-T x - 1.966 x E(B - V) (2) 

We use a foreground reddening value of E(B - V) = 0.11 
for NGC 5128, corresponding to E(C-Ti) = 0.22 for the 
transformation. The uncertainties for a typical color are 
±0.07 dex in the metal-rich and ±0.2 dex in the metal- 
poor regimes. If there is no C-Ti color for the GC, we 
ex amine their U-V c olor and followed the transformation 
of IReed et all (U99l 



[Fe/H]u-i 



-3.061 



(2.015 x (([/ - V) - E(U - V))) 

with a foreground reddening value of E(U — V) — 
E{U - B) + E(B - V) = 0.2. The GCs were 
classified as metal-poor (\Fe/H]c -T, < ~l-0) or as 
metal-rich ([Fe/H] c y, > -1.0) ((Harris et al.l l2004bt 
IWoodlev et al1l2005L 12001) . When neither C, T x , U, 
or V information is available, we divided the GCs as 
metal-rich if (B — I) > 2.072 and as me tal-poor if 
(B - I) < 2.072 following IPeng et al.l (|2004b| ). Our final 
GC sample of 605 has 291 metal-poor, 292 metal-rich, 
and 22 GCs with insufficient photometry for any of the 
above transformations. 

The radial velocity distribution functions are shown 
in Figure [3] for the entire GC population, the red GCs, 
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the blue GCs, and the newly confirmed GCs from this 
study. We fit a Gaussian function to the entire GCS us- 
ing RMI^XQ. We find the best Gaussian fit has a mean of 
517±7kms _1 and a sigma of 160±5 km s -1 . The median 
is v p = 526 km s _1 and the standard Pearson's approxi- 
mation for the mode is v p — 3 x median — 2 x mean = 545 
km s -1 , which is closer to the accepted systemic velocity 
of the galaxy of 541 km s" 1 . The shift to lower mean 
velocity than the galaxy's systemic velocity may suggest 
that our catalog of GCs may still contain contamination 
at the low velocity end from foreground blue Milky Way 
halo stars. We expect to find the higher velocity counter- 
parts to these lower velocity GCs, or hope to remove the 
possible low velocity contamination in future work. We 
do not find any bias in previous radial velocity searches 
which would hinder finding these higher velocity GCs. 

We have analyzed the radial velocity measurements as 
a function of galactocentric distance, shown in Figure |4l 
If a large number of low-velocity objects are Milky Way 
foreground star contamination, we would expect to see a 
large number of blue objects covering this velocity range. 
Looking at Fig. HJ this is not clear, but in the 15'-20' 
range, we can see there are more GCs with lower radial 
velocities than high radial velocities. 

To examine this more closely, we have created veloc- 
ity distributions as a function of galactocentric position, 
shown in Figure [5j We have binned our GCs in units of 
effective radius of the galaxy light (R e // ~ 5'), from 0'- 
5', 5'-10', 10'-15', 15'-20', with the outer GCs in the last 
bin of 20'-45'. We fit the distributions with unimodal 
Gaussians with RMIX, overplotted in the figure, with 
mean, sigma, and reduced x 2 values listed in Table [4] 

While it is clear that the metal-poor GCs on average lie 
at lower radial velocities than the metal-rich GCs, only 
the 10'-15'and 20'-45' bins show the metal-poor GCs lie 
at significantly lower velocities. They may well be af- 
fected by contamination of foreground field stars in these 
outer regions, but the number of known GCs here is quite 
low, and the asymmetries may be purely low number 
statistics. 

3.3. Photometry and Projected Distribution 

With a total sample of 605 GCs in NGC 5128, we 
briefly examine their photometric properties and pro- 
jected distribution. In Figure[6]we show the color, C-Ti, 
as a function of the projected galactocentric radial dis- 
tance in the galaxy. There are 4 very red GCs (GC0078, 
GC0408, GC0411, and GC0552) which by their colors 
may be red background galaxies, but their radial veloc- 
ity measurements are 568 ± 28, 838 ± 63, 548 ± 33, and 
811 ± 32 km s -1 , respectively. The 3 very blue GCs 
(GC0084, GC0282 and GC0432) have radial velocities of 
348±31, 613±6, and 396±84 km s _1 , respectively. These 
blue GCs could either be Milky Way foreground stars or 
they could genuinely be very young GCs. Ages of these 
blue GCs determined spectroscopically are needed to se- 
curely distinguish the difference. We exclude these 7 ob- 
jects from our kinematical analysis below in an attempt 
to examine the kinematics of the normal GC population. 
We have also analyzed the kinematics after removal of 

4 The complete code, available for a variety of platforms, 
is pub licly available from Peter MacDonald' s W eb site at 
http: / /www. math. mcmaster.ca/peter/mix/mix. html 



the 5 GCs that may be stars (discussed in Section l3"7Tj) . 
and find there is no difference of any determined param- 
eter within uncertainties. 

The GC luminosity function is considered a standard 
candle distance indicator. The turn-over magnitude of 
the GC luminosity function of the Milky Way and M31 
GCSs has a mean of V = -7.4 ± 0.2 mag (IHarrisI 11991b 
I Ashman fc Zepl 1x998ft . With our a priori knowledge of 
a distance, we can therefore estimate the turnover mag- 
nitude of the GCS in NGC 5128 to be V = 20.8 mag, 
corresponding to Ti = 20.3 mag. The left hand side 
of Figure [7] shows the luminosity function for 569 GCs 
with the necessary photometry. The distribution appears 
Gaussian with a mean of 19.44 ± 0.04 mag and a sigma 
of 1.03 ± 0.03 mag. Therefore we are clearly missing 
the faintest GCs in the system. We have also fit the 
luminosity functions of the metal-rich and metal-poor 
GC subpopulations and found a mean of 19.28 ± 0.06 
mag and a sigma of 0.90 ± 0.04 mag for the metal-poor 
GCs and a mean of 19.36 ± 0.06 mag and a sigma of 
0.98 ± 0.04 mag for the metal-rich GCs. The metal- 
rich GCs have a slightly fainter peak luminosity than the 
metal-poor GCs, which has also been see n in other GCSs 
(jLarsen et alJ 12001b ITaniura et all 120061 ). Assuming the 
mass distribution between the metal-rich and metal-poor 
GCs is the same (which may not be true if the GC sub- 
populations formed at different times), the difference in 
the turn-over magnitude would be expected if the color 
of a GC reflects its metallicity, which is a good approx- 
imati on for GCs older than a few Gyrs (jWorthev et al.1 
1994). The integrated light of a metal-poor GC would 
appear brighter than a metal-rich GC, at a given age, be- 
caus e the giant stars would b e brighter for a metal-poor 
GC (|Elson fc Santiagdfl^96h 

On the right side of Fig. [71 the color-magnitude dia- 
gram of the GCS shows a clear bimodality, also seen in 
Fig. [BJ Th e photometric uncert ainties for the bright ob- 
jects from iHarris et all (|2004allbT ) range typically from 
0.05 — 0.09 mag, so much of the observed scatter in 
both sequences in Fig. [7j represents real cluster-to-cluster 
metallicity differences. 

The projected radial distribution of the 605 GCs in 
NGC 5128 is shown as a function of projected azimuthal 
angle in the sky (measured in degrees east of north) 
in Figure 1 The GCS of NGC 5128 has been spa- 
tially searched extensively out to 10', while only a few 
searches, including the work presented here, have tried 
to s earch thoroughly ou t to 20'. Beyond this distance, 
only iPeng et al.l (|2004d ) has made an attempt to exten- 
sively search along the isophotal m ajor axis of the galaxy 
of 35° E of N (|Dufour et alJll979l) . Within 20', we do see 
a lack of GCs along the isophotal minor axis (119° E of 
N lDufour et al.lll979[ ). However, with upcoming searches 
for GCs beyond 20', we will be able to test if the GCS 
follows the ellipticity of the underlying galaxy light. 

4. THE KINEMATICS AND DYNAMICS OF THE 
GLOBULAR CLUSTER SYSTEM 

4.1. Mass and Mass-to- Light Determinations 

In the current paradigm of cold dark matter (CDM) 
and hierarchical merging, we expect to find luminous, 
baryonic matter at the centers of much larger dark mat- 
ter halos. Determining the full extent of these dark mat- 
ter halos and the total mass of a galaxy has proven to be 
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quite challenging. 

The technique of measuring integrated stellar light 
spectroscopy is useful for obtaining a velocity dispersion 
from high S/N data in which one can calculate a mass. 
However, the stellar halo gets fainter with galactocentric 
distance, limiting studie s to a few inner effec tive radii 
(jKronawitter et al.ll2000t iGerhard et al.ll200l . HI gas 
has also been successfully used to determine mass con- 
tent, but is inapplicable for early- type galaxies without 
large amounts of gas. Exten ded X-ray halos are found in 
massive early-type gal axies ([Mathews fc Brighentd[2003t 
Humphrey et al.ll2006D which may be used to determine 
the mass out to large distances, where dark matter is ex- 
pected to dominate. This has been d one in luminous e l- 
lipticals such as in NG C 4486 (M87) (|McLaughlinlf l999). 
and NGC 4649 (M60) (|Bridges et al.ll2006t IHwang et al l 

I200I . 

In NGC 5128, tracer particles have been used to 
estimate the mass v ia both PNe (iPeng et al.l 
Woodlev et all 120071) and G Cs (IPeng et al.l 



2004a; 



2004b 



Woodlevl 12006; W oodlev et all 120071) . With our larger 



sample of over 560 GCs with radial velocity measure- 
ments, we can derive an improved estimate of the mass 
and mass-to-light ratio of NGC 5128. 

4.1.1. Pressure and Rotation Supported Mass 

We calculate the mass of NGC 5128 using its GC pop- 
ulation by adding together the component of mass sup- 
ported by pressure, M p , and that supported by rotation, 
M r with 

M tot = M p + M r (4) 

where M p is determined by the Tracer Mass estimator, 
discussed bv lEvans et al.l ((2003) as 



M„ 



(j ^ 



(5) 



where N is the number of objects in the sample, Ri is 
the projected galactocentric radius of the tracer object, 
and Vfc is the radial velocity of the tracer object with 
the rotation component removed. We have determined 
the rotational component for the GCs using the kine- 
matic solution described in Section 14.1.21 The Tracer 
Mass estimator is an advantageous mass estimator over 
the projected and yirial m ass estimators describe d in 
iBahcall fc Tremaind (|198lD and lHeisler et al.l ([1985D be- 
cause it does not assume that the distribution of tracer 
objects follows the underlying dark matter halo. 

We assume the GCs are an isotropic population, so the 
value of C is 



C = 



4(a + 7 )(4 



7)( l_ ( ^)(3-7)) 



(6) 



where Ti n and r out are the 3-dimensional radii corre- 
sponding to the 2-dimensional projected radii of the in- 
nermost, Rim and outermost, R ou t, tracers in the sam- 
ple. The parameter a is set to zero for an isothermal 
halo potential in which the system has a flat rotation 
curve at large distances. The slope of the volume den- 
sity distribution is 7 which is found by determining the 
surface density slope of the sample and deprojecting the 
slope to three-dimensions, written as r -7 . Although the 



Tracer Mass estimator uses a sample of GCs defined be- 
tween ri n and r ou t , it determines the total enclosed mass 
within the outermost point. 

The mass component supported by rotation in Eqn. 2] 
is determined from the rotational component of the Jeans 
Equation, 



M r = 



p „,2 
J^out "max 

G 



(7) 



where R ou t is the outermost tracer projected radius in 
the sample and v max is the rotation amplitude. 

4.1.2. Systemic Velocity, Rotation Amplitude, and Rotation 

Axis 

To obtain the amount of rotation that needs to be re- 
moved from each GC in Equation^ we use the equation 



« p (e) 



J sys 



QRsin(e - G ) 



(8) 

described in lCote et al.l (|200lD . With our known quan- 
tities, the measured radial velocity, v p , and the angular 
position of each GC measured on the projected sky in 
degrees East of North, 0, we fit our GC dataset with 
a weighted least squares fit with the non-linear func- 
tion. From the fit, we extract the systemic velocity, v sys , 
the rotation amplitude (a product of the projected an- 
gular velocity and the projected galactocentric radius of 
the GC) ViR, and the projected rotation axis, O , also 
measured in degrees East of North in the projected sky. 
Eqn. [5] assumes spherical symmetry, that Q is only a 
function of the projected radius, and that the rotation 
axis lies in the plane of the sky. 

4.1.3. Surface Density Profiles 

In order to calculate Equation [51 we require the sur- 
face density profile of the GCS. GCSs typically have ra- 
dial surface density profiles that fall off with a power 
law, with an exponen t that varies between 1.0 — 2.5 
(Ashman fc Zepll 11998). To calculate the surface den- 
sity slope of the GCS in NGC 5128, shown in Figure El 
we bin the GCs into circular annuli with equal number 
of GCs in each to allow near equal statistical weighting 
([Maiz Apellaniz fc Ubedal l2005h . We then fit the data 
with a power-law within the region of approximate az- 
imuthal completeness between 5'-20' for the entire GCS, 
the metal-rich and metal-poor subpopulations of GCs. 
Within 5', a flattening of the GCS is evident in Fig. [9] 
due to the lack of GC candidates in the dustlane of the 
galaxy. Also evident is a drop-off in the GCS beyond 20' 
where there are few to no known GCs along the isopho- 
tal minor axis. Our best fits within these bounds yield 
a volume density slope of 3.38 ± 0.17, 3.26 ± 0.29, and 
3.56 ± 0.21 for the GCS, the metal-poor GCs, and the 
metal-rich GCs, respectively, which are used in Equa- 
tion [6] for all subsequent mass determinations. We also 
determine the volume density profile of the PNe within 
the 5'-20' b ounds to yield a , slope 3.47 ± 0.12 (obtained 
similarly in iWoodlev et al.l 120071 who determined the 
surface density profile excluding only PNe within 5'). 

4.1.4. Mass and Mass-to-Light Results 

Foll owing a similar technique to iSchuberth et all 
(2009), we have used the Tracer Mass Estimator as a 
way to remove GCs that have extreme radial velocities 
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at a large projected distances which would artificially in- 
flate the total mass estimate. This can be seen by their 
contribution to the v 2 R term in Equation First, we 
determine the total mass of the NGC 5128 using all of 
the GCs in our sample. Then, we remove the GC with 
the highest v 2 R contribution from our sample and recal- 
culate the total mass. We continue to calculate the total 
mass while removing the subsequent largest v 2 R contrib- 
utor in each iteration. For each step, we determine the 
change in total mass from the total mass determined us- 
ing the entire GCS, divided by the total number of GCs 
removed. Our results are shown in Figure 1101 We have 
removed the 8 largest v 2 R contributors from the total 
mass estimate, as well as in the kinematic analysis that 
follows. These 8 GCs, 4 of which are metal-rich, are all 
found at a galactocentric radius greater than 16', and 
lead to an overestimate in the mass of NGC 5128. 

To generate a mass profile within NGC 5128 using its 
GCS, we determine the rotation and pressure supported 
masses in cumulative bins from 0'-5', 0'-10', 0'-15', 0'-20', 
0'-45'. The results are tabulated in Table[5]with columns 
displaying the GC subgroup, the radial range of included 
GCs, the number of GCs per radial bin, the systemic ve- 
locity, the rotation amplitude, the rotation axis, the pres- 
sure supported mass, the rotation supported mass, and 
the last column indicates the total mass of the galaxy. 
The total mass determined is the mass enclosed within 
the outermost tracer point. 

We also calculate the mass of NGC 5128 by exclud- 
ing the GCs located within 0-5' projected galactocentric 
radius listed in Table [SJ The Tracer Mass estimator as- 
sumes spherical symmetry of the system and the GCs 
within 5' have severe spatial bias within the dustlane. 
While we examine the enclosed mass out to 45', our most 
secure mass estimate is out to 20' (determined with the 
GCs between 5'-20') to be (5.9 ±2.0) x 10 11 M Q . We also 
determine the total mass using only the metal-poor and 
metal- rich GCs out to 20' and out to 45', also listed in 
Table [5] Both metallicity subpopulations provide good 
agreement for the mass estimate to the entire GCS. 

The mass-to-light ratio of NGC 5128 was also calcu- 
lated using the GCS. We have calculated the absolute 
B-band magnitude assuming a distance of 3.8 Mpc and 
an apparent B-band magnitu de of 7.84 and a galactic ex - 
tinction of 0.5 for NGC 5128 (jKarachentsev et al.ll200l ). 
We determine the B-band mass-to-light ratio, M/Lb, out 
to 20' and to 45' to be 16 M Q /L B o and 30 M /L B o, re- 
spectively. 

4.1.5. Comparison to Previous Measurements 

In NGC 5128, the HI gas shells were used to es- 
timate a mass of 2 x 10 11 M Q within 15 kpc under 
their assumption that NGC 51 28 is at a distance of 3.5 
Mpc (jSchiminovich et al.l 119941 ). With the distance of 
3.8 Mpc used in t h is stu dy, the mass determined by 
iSchiminovich et aH (|1994T> would increase to just less 
than 2.2 x 10 11 M Q within 15 kpc (or 13.8'). Within 
13' we re-evaluate our total mass estimate and find 
M t = (3.5 ± 1.2) x 10 11 M , in agreement with the HI 
g as shells estimate. 

iPeng et al.l (|2004bl ) have used 215 GCs extending out 
to 40 kpc (~ 36') to estim ate a pressure suppo rted 
mass of 7.5 x 10 11 M (see IWoodlev et ail [20071 for 



this mass, corrected from IPeng et al.l ([2004b)). We re- 
evaulate the pressure supported mass out to 36' and ob- 
ta in M„ = (9 7 ±3.3) x 10 11 M , matching the estimate 
o f IPeng et al.l (l2004bl) w ithin uncertainties. 

IWoodlev et al.l (|2007| ) found a total mass, including the 
contribution from rotation and pressure, of (1.3 ± 0.5) x 
10 12 M Q from 340 GCs out to the full extent of the known 
GCS (45'). We find M t = (1.1 ±0.4) x 10 12 M Q in excel- 
lent agreement with this result although we use a much 
larger s ample of 429 GCs (e xcludi ng GCs with R„ c < 5') . 

Both IPeng et al.l (|2004bf ) and IWoodlev et ail ((20071 ) 
have used samples with quite low number of GCs be- 
yond 10' in galactocentric radius. We do, however, find 
good agreement with their earlier work, but with bet- 
ter constraints on the uncertainties associated with the 
determined mass. 

4.2. Kinematic Analysis 

The previous kinematic studies of the GCS in NGC 
5128 have been done with smaller samples of GCs, with 
215 GCs used in lPeng et "all (l2004bl\ and 34 GCs used in 
IWoodlev! (f2006h and IWoodlev et al.l (|200l . Here, with a 
larger sample size and better spatial coverage, we evalu- 
ate the kinematics of the GCS. 

We determine the kinematics for the GCs using Equa- 
tion [8] for the full sample or discretely binned subsamples 
by weighting the individual GCs according to their indi- 
vidual radial velocity uncertainties, v p ^ err and the ran- 
dom velocity component of the system, v ran . The ra- 
dial velocities and associated uncertainties used in this 
study are the weighted averages of all previous measure- 
ments. The random velocity component is the standard 
deviation of the radial velocities of the GCs from the 
best fit rotation curves determined without any weight- 
ing. The final uncertainty for each GC used in Eqn. [5] 
is a = (v 2 err + v 2 an y/ 2 , and the weight is 1/er 2 . The 
dominant uncertainty is v ran providing all GCs a near 
equal weighting in the fitting. 

Our initial results are shown in Figure [TT] for the entire 
GCS, and for the metal-poor and metal-rich subpopula- 
tions of GCs. The best sine curve fits clearly show the 
large dispersion (or random motion) of the GCS with a 
small rotational component. Our results are presented in 
Table [6] which lists the GC subgroup, the radial range of 
the GC coverage, the average radius, and number of GCs 
within the subgroup, followed by the best fit systemic 
velocity, rotation amplitude, r otation axis, the is opho- 
tal major axis of 215°E of N (|Dufour et al.lll979f ) sub- 
tracted from the rotation axis, velocity dispersion, and 
rotation parameter (rotation amplitude/velocity disper- 
sion). The latter two will be discussed in Section 14.2.11 
We note the kinematics presented in Table [5l while cal- 
culated with the same method as in Table El differ by 
their radial ranges. In the mass determination presented 
in Section 14. 1[ we were developing a mass profile that 
with the inclusion of more datapoints reduces the uncer- 
tainties on the mass determined using the enclosed data. 
Here, we are searching for the kinematic signatures at 
distinct radial positions. 

For the entire GCs, our best fits are v sys = 517 ± 7 
km s~\ OR = 33 ± 10 km s" 1 , and o = 185 ± 15°E 
of N. Also, our best fits for the overall kinematics are 
v sys = 506 ± 9 km s"\ VIR = 26 ± 15 km s" 1 , and O = 
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177 ± 28°E of N for the metal-poor and v sys = 526 ± 10 
km s" 1 , OR = 43 ± 15 km s~\ and 6 D = 196 ± 17°E of 
N for the metal-rich GC subpopulations. As discussed 
in Section 13.21 the overall velocity of the metal-poor GC 
subpopulation is shifted to a lower value than the metal- 
rich GC subpopulation, however, both are lower than the 
systemic velocity of the galaxy. The rotation amplitude 
as first indicated in Fig. [Tl] is quite low. The metal-poor 
GC subpopulation has very mild rotation and the metal- 
rich GCs have mild to moderate rotation as a whole. The 
metal-rich GC subpopulation appears to rotate around 
the isophotal major axis of the galaxy within 15', while 
the metal-poor GC subpopulation does not appear to 
follow either the isophotal minor or major axis. We note 
there is no difference in the results within uncertainties 
if v sys is fixed at the systemic velocity of the galaxy or 
left free in Eqn. [8] 

We create kinematic profiles for the entire GCS as well 
as the metal-poor and metal-rich GC subpopulations as 
functions of galactocentric radius. We bin the GCs in 
galactic effective radial bins of 0'-5', 5'-10', 10'-15', 15'- 
20', and the remaining GCs in the final bin of 20'-45' for 
the entire GCS and the metal-poor GC subpopulation. 
The metal-rich GC subpopulation has a more spatially 
centrally concentrated system so we bin in 0'-5', 5'-10', 
10'-15' bins, with the remaining metal-rich GCs in the 
final bin of 15'-45'. 

We list the results in Table [6] for the best fit systemic 
velocity, rotation amplitude, and rotation axis for the 
entire GCS, as well as metal-poor and metal-rich GC 
subpopulations. With radial distance, we see both the 
metal-rich and metal-poor populations maintain a rela- 
tively constant amount of rotation with radial distribu- 
tion, with the metal-rich GCs having a larger rotation 
amplitude than the metal-poor GCs. Note that the ro- 
tation amplitude for the entire metal-poor GC subpop- 
ulation from 0'-45' is less than most of the values deter- 
mined in the specific radial bins,which at first sight is an 
anomaly. However, this curious result is connected with 
the solutionos for the rotation axis angles, which differ by 
large amounts from bin to bin. When the rotation am- 
plitude is zero, the rotation axis is ofcourse undefined, so 
when we average together the results from the different 
radial bins (all with a different axis of rotation), they 
average down to zero with no meaningful axis. We con- 
clude from this that the rotation in the metal-poor GC 
subpopulation is very mild. There is an indication that 
the metal-poor GCs have an increased amount of rotation 
in the outer regions (beyond 20'), but again, the GC pop- 
ulation in this region clearly violates the basic assump- 
tion of spherical symmetry used in Eqn. [8] We therefore 
caution as to its validity. If the metal-poor GC system 
rotation does increase, this may support the formation 
mechanisms of disk-disk mergers for NGC 5128. Here, 
the orbital angular momenta of the progenitor galaxies 
are converted to intrinsic angul ar momentum in t he rem- 
nant from the merging event (jBekki et al.l 120051 ) . This 
angular momentum is transferred to the outer regions of 
the galaxy, and thus it would be expected for the rotation 
amplitude of the GCS to increase at larger radii. 

The metal-poor GC subpopulation shows a potential 
twist in its rotation axis within the inner 5' from ro- 
tation around the isophotal minor axis to around the 



isophotal major axis. Beyond 10' the axis may change 
again. However, given the near zero rotation amplitude 
of the metal-poor GCs, it does not appear that they ro- 
tate around a well-defined rotation axis. In the inner 20' 
the metal-rich GCs are consistent with rotation around 
the isophotal major axis of the galaxy. In the outer re- 
gions, the rotation axis of the metal-rich GC subpopu- 
lation appears to decrease slightly, but again, this is a 
region of uncertainty. 

4.2.1. Velocity Dispersion 

The velocity dispersion of the system is determined us- 
in g a maximum like l ihood dispersion estimator discussed 
in lPrvor & Mevlanl (|1993| ). 



N 
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where N is the number of GCs in the sample, Vi is the 
radial velocity of the GC after subtraction of the rota- 
tional component determined in Eqn. [8] with v sys held 



constant at 541 km s 1 , and 



is the uncertainty in 



the velocity measurement. The projected velocity dis- 
persion, a vp, is determined by iteration until satisfying 
Equation [3] The systemic velocity used here was deter- 
mined from Eqn. |SJ The uncertainties in the velocity 
dispersion are c alculated from the vari ance of the dis- 
persion, l isted inlPrvor fc Mevlanl ([19931 ) . This has been 
shown bv lPrvor fe Mevl an (1993) to estimate the uncer- 
tainty to 5% with as little as 20 objects in the sample. 

We list the velocity dispersions in Table [51 As a 
whole the metal-poor and metal-rich GC subpopulations 
have identical velocity dispersions within uncertainties 
of 149 ± 4 km s" 1 and 156 ± 4 km s™ 1 , respectively. 
However, the behaviour with galactocentric radius may 
be different for the different metallicity subpopulations. 
The metal-poor GC velocity dispersion is flat or declining 
out to 20' beyond which it appears to rise. The metal- 
rich GC profile shows a declining velocity dispersion out 
to 5' after which it appears to flatten, which may indi- 
cate the metal-rich GCs are more radially biased than 
the metal-poor GCs. A rise in the velocity dispersion 
in the outer regions of a galaxy could be an indicator 
for large amounts of dark matter or of anisotropic orbits 
of the velocity ellipsoid of the GCS itself. With an ex- 
tended X-ray gas profile to trace the mass out to regions 
beyond the halo light, we could determine the degree of 
anisotropy of these orbits, if they ex ist, using a Jean s 
soluti on. This has been done fo r M87 (iCote et alll2001l) . 
M49 dCote et al.l 120031). M60 dHwang et al.l 120081). and 
NGC 1407 (|Romanowksv et all I2009D as examples. It 
may be possible to find indications of anisotropy in the 
GCS by examining its radial velocity as a function of ra- 
dius. To analyse the ellipsoid distributions, it is more 
instructive to remove the rotational component of the 
GC subpopulations; however, considering the rotation 
amplitudes are nearly negligible, we use Fig. [5] as a first 
glance into possible anisotropy of the GCS. 

If the GCs are moving on very radial orbits, we expect 
the distribution to be strongly peaked at the systemic 
velocity, while if they are moving on very circular or- 
bits, the distribution would have a more flattened shape. 
Isotropic orbital motion would have a Gaussian distribu- 
tion. The radial velocity distributions do not appear to 
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be tightly peaked in any histogram in Fig. [3J but there 
is a slight tendency for the GC velocity histograms to 
have a wider distribution at larger galactocentric radius, 
thus a flattening effect. In Tabled] we quote the reduced 
X 2 values for these Gaussian fits. We note a general 
trend for larger Xred values for the outer bins, indicat- 
ing a poorer fit of the model to the data, perhaps by 
the lower number of GCs in the outermost regions. A 
flattened distribution may also be caused by the con- 
tribution of GCs from accr eted satellite gala xies in the 
outer regions of the galaxy (jBekki et al.l l2003). Evidence 
for potential substr ucture has been presented in M31 
(IMerrett et al.1 12001 . M87 (IDohertv et alJ l2009h. NGC 
1399 (ISchuberth et alJl2009h . and perhaps in NGC 1407 
(|Romanowksv et al.ll2009f )~Interestinglv. numerical sim- 
ulations have suggested that tracer populations should 
have isotropic orbits in the inner regions of their host 
galaxy, while in the outer regions, the i r orbits should 
become anisotropic (iDekel et al.l 120051: iDiemand et al.l 
[20051: lAbadi et al.l 120061: IGnedin fc Prietol 120061) . In re- 
lation to the GC results presented here, the presence of 
substructure of accreted satellites or interlopers will bias 
the velocity dispersion to higher values. This is seen 
in the outer-most bin in the metal-poor GC subpopula- 
tions, possibly indicating substructure from an accreted 
satellite or Milky Way halo star contamination. 

The rotation parameter is a measure of the maxi- 
mum rotation divided by the velocity dispersion in a 
system. We obtain rotational parameters of 0.22 ± 0.07, 
0.17 ± 0.09, and 0.28 ± 0.10 for the entire GCS, metal- 
poor, and metal-rich subpopulations of GCs, respectively 
(listed in Table [S]). These values reinforce what we have 
already discovered about the GCS. While the velocity 
dispersion of the entire metal-poor and entire metal-rich 
GC subpopulations is consistent within uncertainties, the 
metal-rich GCS, as a whole, has more rotation. 

4.2.2. Kinematics for Globular Clusters with Ages 

iWoodlev et al.l (|2010f) have estimated ages, metallici- 
ties, and a-to-Fe abundance ratios for 72 GCs in NGC 
5128. By calibrating their high signal-to- noise spec- 
troscopy tak en with Gemini-S/GMOS to the Lick in- 
dex system dBurstein et al.| H9841: IWorthev et al.l 119941: 
IWorthev fc Ottavianilll997l: iTrager et al.lll998D .lhev ex- 
tracted their estimates by comparing their measured 
indices to th e simple st ellar population models of 
iThomas et all (|2003l I2004Q . Using Eqn. we have es- 
timated the kinematic properties of these 72 GCs, as 
well as subsets of this group, such as their GCs esti- 
mated to be old (ages > 8 Gyr), intermediate (5 — 8 
Gyr), and young (age < 5 Gyr). We also analyse the 
old GCs in subgroups of metal-rich and metal-poor. Our 
results are presented in Table [7l which lists our subgroup 
of GCs, the number of GCs in each subgroup, the sys- 
temic velocity, the rotation amplitude, and the rotation 
axis. The kinematics for the GCs subdivided by metallic- 
ity and age have quite large uncertainties because they 
are small samples. They do, however, indicate general 
properties of these small subsets of GCs. As expected, 
the old metal-poor GCs and old metal-rich GCs generally 
match the kinematic properties of the entire metal-poor 
GCs and metal-rich GCs, respectively. The rotation axis 
of the young GCs, which are all metal-rich, is quite dif- 
ferent than that of the entire metal-rich GC population. 



With a rotation axis of 80 ± 84° E of N at the average 
projected radius of 5.7' this is significantly different from 
the metal-rich GC rotation axis of 244 ±27° E of N at an 
average radius of 3.49' or 178 ± 22° E of N at an average 
7.23' distance. The young metal-rich GCs appear to be 
rotating either along the isophotal major axis, but in the 
opposite direction to the metal-rich GC subpopulation, 
or along the isophotal minor axis, or anywhere in be- 
tween. In any case, their rotation axis is quite different 
from the bulk of the metal-rich GCs. 

4.2.3. Kinematics for Globular Clusters versus Magnitude 

The kinematics of the GCS with respect to Ti magni- 
tude has also been analyzed upon removal of GCs that 
have extreme colors (see Section I3.3[) . The GCs are 
binned in magnitudes between 16—18, 18—18.5, 18.5—19, 
19.5 — 20, 20.5 — 23 mags, in order to keep the number 
of objects per bin above a minimum of 25. We plot the 
results for the systemic velocity, rotation amplitude, ro- 
tation axis, and velocity dispersion in Figure [T2l We find 
there are no major differences in the kinematic properties 
of the GCs with changing magnitude. It does appear, 
however, that the GCs between 18 — 18.5 magnitudes 
have a significantly different rotation axis (41 ± 34° E of 
N) than the remaining GCs which have an axis of ~ 200° 
E of N. We have examined the projected distribution of 
the GCs in this bin compared to the other magnitude 
bins and see no distinct difference. We also notice that 
the brightest GCs have a significantly larger velocity dis- 
persion than the remaining GCs, indicating there may 
be evidence for varying kinematics as a function of mag- 
nitude. A similar finding was identified in NGC 1407 
(|Romanowksv et all 120091 ): however, as they indicated, 
it could be due to spatial bias of the brighter objects in 
the sample being close to the center of the galaxy. This 
is not as obvious in the case of NGC 5128, as almost all 
GCs in this study are centrally concentrated within 15'. 

4.2.4. Kinematics of the Planetary Nebulae 

The PNe are a field star population at one of the lat- 
est stages of stellar evolution and are easily identified by 
their emission of [OIII] . Follow-up spectroscopy can then 
be performed on the PNe to measure their radial veloc- 
ities which has been done for 780 PNe in NGC 5128 
extend ing out to over 80', and compiled in iPeng et all 
(|2004aft . 

The k inematic prope r ties o f the 780 PNe are investi- 
gated in lWoodlev et al.l (|2007l ) in the same manner as the 
GCs presented in this study . They found the PNe had 
a rotation amplitude of 76 ± 6 km s , a rotation axis of 
170 ± 5° E of N, and a velocity dispersion of 118 ± 13 km 
s _1 . The rotation amplitude of the PNe is a bit larger 
than the GCS as a whole and the rotation axis of the 
P Ne is quite uniform w ith increasing radius (see Table 5 
in lWoodlev et al.l l2007). The rotation axis is also, inter- 
estingly, similar to the isophotal major axis of the galaxy, 
closely mimicing the metal-rich GC subpopulation in the 
inner 15'. The velocity dispersion of the PNe is different 
than the GCs. The PNe are dynamically colder in the 
inner 5'and their velocity dispersion continues to decline 
with increasing galactocentric radius. 

Declining velocity dispersions of the PNe with 
galact ocentric radius have b e en found for NGC 
3379 (|Romanowskv et all l2003t iDouglas et all 120071) , 
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NGC 821 (IRomanowskv et al.1 l200l. NGC 4494 
(IRomanowskv et al.ll2003t iNapoliano et al.1 12009ft . NGC 



5128 (IPeng et alJl2004aHWoodlev et al.ll2007D . and NGC 
4697 (jMendez et al.l 120091 ) . which could indicate a lack 
of dark matter in the system (jRomanowskv et al.l[2 003). 
Howe ver, GCSs studied within NGC 3379 (jPierce et al.1 
2006) have shown an increasing velocity dispersion, indi- 
cating there is likely to be some amount of dark matter 
present . Rec ent work has shown that PNe may be on 
radial orbits (iDekel et al.l 120051 : iMamon fe Lokasl 120051 : 
lAbadi et all 12006). which may also cause the declining 
velocity dispersions. Anisotropy may not necessarily 
be expected of the GCs because these two populations 
trace different kinematical histories of their host galax- 
ies, especially if the galaxy has undergone a major merg- 
ing event in the past. Other studies have found that 
the GC populations have approximately isotropic orbits 
(iZepf et al.l [20001: ICote etall [2001 l200l Bridges et afl 
120061: iSchuberth et al.l 12009ft . However, separating the 
blue and red GCs have s hown there may be slight 
anisot ropy as well, with M87 (ICote et al.l|200lT) and NGC 
4649 (|Bridges et al.l l2006t iHwang et alM OOsI as exam- 
ples. In the case of NGC 5128 studied here, we therefore 
suggest that the PNe may be following more radial orbits 
than the GCs. 

5. DISCUSSION 

The kinematics of the GCSs that have been performed 
in other giant elliptical galaxies base d on current samples 
of > 100 GCs includes: NGC 1407 (jRomanowksv et al.1 
2009), an elliptical gala xy in the Eridanus A group at 
a dista nce of 20.9 Mpc dForbes et al.|[2006h: N GC 4649 
(M60) (jBridees et al.l 120061: IHwang et al.l 12001 . a lumi- 
nous giant elliptical galaxy i n the Virgo clust er of galaxies 
at a dis tance of 17.3 Mpc dMei et all 12007ft: NGC 4486 
CM871 dKissler-Patig fc Gebhardtl 119981: ICohenl l2000l : 
ICote et al.H200l 7 a superm assive cD galaxy , located at a 
distance of 16.7 ± 0.2 Mpc (jMei et alj | 2007h in the Virgo 
clust e r of galaxies; N GC 4472 (M49) ([Sharpies et al.1 
[19981 IZepf et al.l l2000l : ICote et all 120031 ) . a giant ellip- 
tical galaxy and also the brightest member of the Virgo 
clust er of galaxies, l ocated at a distance of 16.4 ± 0.2 
Mpc (IMei et alObol : NGC 4636 (ISchuberth et alfEM 
iLee et al.ll2010ft . the southernmost giant elliptical galaxy 
from the dynamical center of the Virgo clust er of galax- 
ies, located 14.7 Mpc away (jTonrv et alJl2001ft : and NGC 
1399, a giant elliptical galaxy and also the brightest 
galaxy in the Fornax cluster of galaxies, 1 9 Mpc away 
(jRichtler et al|[200l ISchuberth et al.l 12009ft . The main 
findings of these st udies have been well summariz ed in 
IHwang et all (|2008ft and iRomanowksv et al.l (|2009ft and 
we do not repeat these here. 

We find that the kin ematic trends of NGC 1407 
(jRomanowksv et al.ll2009l ) are quite similar to NGC 5128 
found in this study. Both systems have larger rotation 
in the metal-rich GCs (at least in the inner-most re- 
gion where both datasets have strong azimuthal cover- 
age), with the rotation of the metal-rich GCs around 
the isophotal major axis. Both systems have velocity 
dispersions that appear flat or decreasing in the inner 
regions with a potential increase at larger radii. They 
are also both dominated by dispersion over their entire 
GC population. Both NGC 5128 and NGC 1407 are 
the only two elliptical galaxies in a group environment 



with large samples of GCs available in order to perform 
a kinematic analysis. With a sample of 2 galaxies, it is a 
far stretch to propose that the GCSs in elliptical galax- 
ies found in group environments have similar kinematics 
due to similar formation histories. With this in mind, 
if these galaxies hierarchically built up the majority of 
their stellar content at early times, later small accretions 
and minor merging will not severely disrupt the orbital 
and kinematic properties in the outer halo. The kine- 
matics from the major episodes of star formation within 
these galaxies may still be intact. 

5.1. Implications for Galaxy Formation 

Th e monol i thic collapse scen a rio (lEggen et all 
19621 iTinslevI [1971 ILarsonl Il97l H975I: ISilld 119771 



Arimoto fe Yoshiil 11987ft suggests the GCs in massive 
galaxies should be old, forming at high redshift. The 
collapse would be very rapid indicating that the spread 
in old ages should be small and the kinematic properties 
of the GCs should be quite similar. Any rotation that 
exists in the GCs would come from tid al torques from 
nearby forming galaxies (jPeebles! 119691 ) . The expected 
amount of rotation would be small, as the angular mo- 
mentum of the initial cloud would be small to begin with 
in order to form the spheroidal shape of the elliptical 
galaxy. This scenario is not supported by the trend 
towards younger ag e s as metallicity increases in the GCS 
(IPeng et all I2004ri IBeaslev et all 120081: IWoodlev et"aH 
12010ft . Neither is it supported by the large amounts of 
rotation that are seen for some GCSs (for the metal-rich 
GCs in NGC 4636, and both metal-rich and metal-poor 
GCs in NGC 4486 and NGC 4649) or by the different 
rotational properties of some metal-rich and metal-poor 
GCs (in NGC 4636, NGC 4486, and NGC 4472, for 
example) . 

The galaxy mer ging model (jSchweizerl 119871 : 
lAshman fc Zepj["l992f) would naturally suggest younger 
ages of the metal-rich GCs, which form, according to 
this scenario, only in the merging events, while the 
metal-poor GCs pre-exist in the progenitor galaxies. 
This theory expects a multimodal c olor distribu t ion o f 
the GCS. Numerical simulations by Bckki ct al. (2002) 
have shown that the metal-poor GCs in this scenario are 
spatially extended while the newly formed metal-rich 
GCs are more centrally concentrated and more extended 
along the isophotal major axis of the galaxy. While 
the above properties are not disputed here, this does 
not explain the old ages for the majority of metal-rich 
GCs found in previous work on NGC 5128, unless all 
the major merging events took place at very early times 
(~ 10 Gyr ago). Any existing angular momentum 
would be transferred to the outer reg ions of the remnan t 
galaxy from the merging process ([Bekki et al.l [2005). 
The more extended metal-poor GC population should 
have a higher rotation signature than the metal-rich 
GCs, particularly in the outer regions. We do not see 
this in the GCs of NGC 5128 or in the GCs of the other 
giant elliptical galaxies studied, with the exception 
of NGC 4472, which generally agrees with the galaxy 
merger model. 

A hierarchical merging scenario, co mbining the 
multi phase dissipational collap se model (Forbes et all 
1997ft and accretion scenarios (ICote et al.1 119981 120001 
2002|i . proposed by IBeaslev et al.l (|2002i 120031 ) and 
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IStrader et al.l (|2005l ). is more consistent with our find- 
ings for NGC 5128. This naturally explains the GC color 
bimodality and also predicts old ages for both metal- 
rich and metal-poor GCs (with the metal-poor GCs hav- 
ing slightly older ages). The old ages obtained for the 
stellar population in the halo of NGC 5128, as well as 
similar chemical enrichment between the metal-ric h GCs 
and stellar halo (|Reikuba et al J 120051: Harris! I2008D . are 
consistent with this scenario as well. This hierarchical 
scenario also supports the younger metal-rich GCs in 
NGC 5128 forming in either more recent merging events 
and/or on a more extended timescale. In the sparse 
environment of NGC 5128, the formation of the more 
metal-rich GCs could be extended because the density in 
the group environment is lower than that in clusters of 
galaxies. This is supported by both the range of ages for 
the metal-rich GCs as well as the [a/Fe] values suggest - 
ing longer formation timescales (jWoodlev et al.1 12010). 
compared to GCs in other giant galaxies. In a typical 
massive collapse, the metal-poor GCs would not neces- 
sarily be expected to form with large amounts of rota- 
tion and instead they sh ould have a large velocity disper- 
sion (|Hwang et al.ll2008l ). However, with the accretion of 
small protogalaxies in the early Universe, the GCs may 
have some rotation as well as orbital anisotropy. The 
metal-rich GCs, in this formation scenario, can vary in 
rotation, depending on the degree of the second collapse 
(jForbes et al.l [1997). We see evidence for this in NGC 
5 128 in both t h e old ages and kinematic results. 

iBekki et al.1 (|2005h have numerically simulated the 
kinematic signatures of GCs from the merging of Milky 
Way-type disk galaxies with pre-existing old metal-poor 
and metal-rich GCs. To simplify their simulations, they 
assumed that the initial GCSs in the progenitor galaxies 
were supported by dispersion with no ro tation. How- 
ever, we see rotation in the GCSs of M3 1 (|Perrett et al.l 
[2001 and the Milky Way (Harris! I200ll ). They also as- 
sumed the mergers were dissipationless, but we do see 
some small fraction of young GCs in NGC 5128. How- 
ever, these numerical simulations allow us to compare 
the kinematic properties that we see in the GCS of NGC 
5128 to the generally expected kinematics from two old 
merging spiral galaxies. Their results indicate that both 
the metal-rich and metal-poor populations would have 
incre asing rotational si gnatures with galactocentric ra- 
dius. IBekki et al.l (|2005l ) also find the velocity dispersions 
of the GCSs are flat or declining out to large radii caused 
by their radial anisotropic orbits. Lastly, they found a 
variety of kinematic alignments for the metal-rich and 
metal-poor GCs. 

In our kinematic study of NGC 5128, we do not see 
a strong case for significant rotation in either metallic- 
ity subpopulation, and also no significant increase with 
radius in the rotation signature. We do, however, see 
indications for flat and declining velocity dispersions in 
the inner 20' of NGC 5128 for both metal-rich and metal- 
poor GCs, and a different axis of rotation for the metal- 
rich and metal-poor GCs. The results found here are 
not clearly consistent with the idea of the GCs forming 
their kinematic signature from a major disk-disk merger, 
though this kind of origin is not ruled out. In the case of 
other galaxies, NGC 1399 shows little to no rotation in 
either metal-rich or metal-poor GCs and NGC 1407 also 
show little to moderate rotation signatures as well. 



There is an impo rtant implication from the study 
of IBekki et al.l (|2005| ) that should not be underplayed. 
Their results suggest that the initial kinematic signa- 
tures of GCs in the progenitor galaxies undergo orbital 
mixing and it may, therefore, not be possible to trac e 
the original kinematics of the GCs (jHwang et al.l f2008'). 
It may only be possible to examine the GC kinematics 
from the most recent ma j or interaction. As suggested by 
iKissler-Patig et all (|1998l) , it would then be very difficult 
to use the current kinematic signature to trace the orbital 
history of GCs. In this respect, we can suggest that the 
most recent major interaction in NGC 5128 was not likely 
a major disk-disk merger. Again, it seems more likely 
that the galaxy formed hierarchically in early times, with 
many inf ailing protogalaxies. 

6. CONCLUSIONS 

We have obtained spectroscopy with LDSS- 
2/Magellan, VIMOS/VLT, and Hydra/CTIO in order 
to measure the radial velocities of GCs in the nearby 
elliptical galaxy, NGC 5128. With these datasets, we 
have remeasured radial velocities for 218 known GCs 
as well as confirmed 155 new GCs. There are now 605 
GCs confirmed in NGC 5128 and 564 of these GCs have 
measured radial velocities. This is the second largest 
kinematic dataset of G Cs in any galaxy after NGC 1399 
(|Schuberth et al.ll2009h . 

The large sample of GCs with radial velocity measure- 
ments allows us to perform a new kinematical analysis of 
the system. We have investigated the systemic velocity, 
the rotation amplitude, the rotation axis, and the veloc- 
ity dispersion for the entire GCS, as well as the metal-rich 
and metal-poor GC subpopulations. We have also ana- 
lyzed these kinematic properties as a function of galac- 
tocentric radius. We find the metal-poor GCs have a low 
rotation 26 ± 15 km s -1 extending out to 45', and does 
not have a well-defined rotation axis. The motion of the 
metal-poor GCs is supported by dispersion with a mea- 
sured value of 149 ±4 km s -1 . It has a flat/declining ve- 
locity dispersion profile extending to 20', at which point 
it appears to increase. The metal-rich GCs have small to 
moderate rotation of 43±15 km s _1 out to 40' around the 
isophotal major axis of the galaxy. It has similar velocity 
dispersion for the whole system (156 ± 4 km s _1 ) com- 
pared to the metal-poor GCs. The metal-rich velocity 
dispersion profile decreases out to a galactocentric ra- 
dius of 5' and then flattens. This increase in the velocity 
dispersion profiles in the metal-poor GC subpopulation 
may be a result of the incomplete and spatially biased 
sample beyond 15' in the GCS, or caused by Milky Way 
halo contamination. Or perhaps it may be due to the 
presence of a substructure among the GCS which could 
artificially inflat the velocity dispersion in the outer halo, 
or it could also be an indicator of an extensive dark mat- 
ter halo. As a consequence of the low rotation amplitude, 
both the metal-poor and metal-rich subpopulations have 
low values of the rotation parameter of 0.17 ± 0.09 and 
0.28 ±0.10, respectively. 

We have also analyzed the kinematics for the GCs that 
have estimated ages and metallicities from the work of 
iWoodlev et all (|2010D and found that the young (< 5 
Gyr) metal-rich GCs rotate on a different axis than the 
entire metal-rich GC population. If these GCs were 
formed from the accretion of a gas-rich satellite galaxy, 
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we could expect their kinematic signature to be different 
from the older GCs. 

We estimate the mass and mass-to-light ratio of NGC 
5128 using its GCs. We find a mass of (5.9 ± 2.0) x 
10 11 Mq out to 20' with M/L B = 16.4 M Q /L BQ . Our 
mass estimates determined using the metal-rich GCs and 
the metal-poor GCs all agree within uncertainties. This 
mass estimate clearly puts NGC 5128 in the class of giant 
elliptical galaxies, with a potentially large dark matter 
halo. 
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TABLE 2 

Radial Velocity Measurements of Previously Confirmed Globular Clusters in NGC 5128 
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395±67 


447±49 




426±21 


GC0144 


13 


25 


11.98 


-43 


04 19.3 


865±77 


612±11 




622±18 


GC0145 


13 


25 


12.11 


-42 


57 25.2 




487±12 




490±10 


GC0146 


13 


25 


12.21 


-43 


16 33.9 


510±25 h 






573±16 


GC0148 


13 


25 


12.45 


-43 


14 07.4 




467±24 




526±33 


GC0149 


13 


25 


12.84 


-42 


56 59.8 




586±29 




664±141 




IS 


9R 


19 QD 


-43 


0.7 RQ 1 




finn4-ia 
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TABLE 3 

Radial Velocity Measurements of Newly Confirmed Globular Clusters in NGC 5128 



GC ID 


R.A. 


Decl. 


C a c 


M 


cm 


Ti 


°Ti V r> LDSS-2 Vr,VIMOS v r, Hydra v r 




(J2000) 


(J2000) 


(mag) (mag) 


(mag) 


(mag) 


(mag) 


(mag) (km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) 



GC0416 


13 


24 


33 


.63 


-43 


12 


01 


6 


21 


,69 





.02 


20 


,94 


0.01 


20 


.12 


0.01 


GC0417 


13 


24 


42 


.39 


-42 


51 


49. 


1 


20 


,50 





,02 


19 


.82 


0.02 


19 


.09 


0.01 


GC0418 


13 


24 


41. 


.77 


-43 


06 


33. 


3 


22 


,01 


0, 


,06 


20 


.86 


0.05 


19 


,87 


0.03 


GC0419 


13 


24 


55 


.31 


-43 


10 


39. 


3 


21 


,87 





,02 


20 


.83 


0.02 


19 


,89 


0.01 


GC0420 


13 


24 


55 


.97 


-43 


02 


15 


9 


21 


.29 





.03 


20 


.64 


0.02 


19 


.87 


0.03 


GC0421 


13 


24 


59. 


.58 


-43 


06 


40. 


8 


22 


.62 





.03 


21 


.49 


0.03 


20 


,47 


0.01 


GC0422 


13 


25 


03 


.28 


-43 


08 


14 


4 


21 


.11 





.03 


20. 


.56 


0.02 


19 


,76 


0.01 


GC0423 


13 


25 


04. 


.28 


-43 


17 


10 


.8 


22 


.36 





.01 


21 


,80 


0.03 


20 


.92 


0.02 


GC0424 


13 


25 


06 


.87 


-43 


02 


10. 


2 


20 


.65 





.01 


19 


,93 


0.01 


19 


.13 


0.01 


GC0425 


13 


25 


09 


.60 


-43 


04 


37 


9 


21 


.65 





,02 


20 


,94 


0.01 


20 


.11 


0.01 


GC0426 


13 


25 


14 


18 


-43 


04 


46 


6 


21 


,23 


0, 


,03 


20 


,67 


0.02 


19 


.87 


0.02 


GC0427 


13 


25 


21 


.49 


-43 


19 


29. 


6 


19. 


.59 


0, 


.01 


19 


.23 


0.01 


18 


.66 


0.01 


GC0428 


13 


25 


22 


16 


-43 


16 


38 


.8 


21 


,06 





.04 


20 


.48 


0.04 


19 


.76 


0.02 


GC0429 


13 


25 


22 


.88 


-43 


08 


10 


1 


20. 


,62 





.01 


19 


.20 


0.01 


18 


.14 


0.01 


GC0430 


13 


25 


24. 


.38 


-42 


57 


16 


.8 


- 




- 




19 


.25 


0.29 


18 


,90 


0.01 


GC0431 


13 


25 


28. 


.72 


-42 


50 


12 


3 


20 


,60 





.04 


19 


.93 


0.04 


19 


,36 


0.01 


GC0432 


13 


25 


34 


.08 


-43 


10 


45 


,7 


19 


,87 





.31 


20 


.81 


0.07 


20 


,00 


0.05 


GC0433 


13 


25 


36 


.11 


-42 


56 


16. 


,8 


22 


.61 





.03 


22 


,58 


0.03 


22 


.00 


0.09 


GC0434 


13 


25 


37 


.85 


-42 


56 


28. 





20 


,20 


0, 


.01 


20 


,59 


0.01 


18 


.51 


0.00 


GC0435 


13 


25 


38. 


.53 


-42 


57 


19. 


9 


21 


,32 


0, 


.02 


20 


.30 


0.03 


19 


.42 


0.01 


GC0436 


13 


25 


41 


.71 


-42 


58 


31 


6 


20. 


,34 





.02 


19 


.29 


0.01 


18 


,52 


0.01 


GC0437 


13 


25 


45 


.23 


-42 


54 


23 


,8 


22 


,04 





.03 


21 


.00 


0.04 


20 


,14 


0.02 


GC0438 


13 


25 


45. 


.91 


-42 


51 


41, 


2 


21, 


.34 





,02 


21 


.55 


0.01 


19 


.56 


0.01 


GC0439 


13 


25 


47. 


.00 


-42 


55 


29 


9 


20. 


,56 





.05 


19 


.75 


0.03 


18 


.89 


0.02 


GC0440 


13 


25 


50. 


.59 


-42 


51 


40 


1 


20 


.19 


0, 


.02 


19 


.60 


0.01 


18 


.95 


0.02 


GC0441 


13 


25 


54 


.48 


-42 


57 


27, 


6 


20 


.97 





.02 


19 


.91 


0.01 


19 


.12 


0.01 


GC0442 


13 


26 


09 


.37 


-42 


53 


17, 


,5 


21 


.09 





.03 


20 


,38 


0.02 


19 


.67 


0.02 


GC0443 


13 


26 


12 


.11 


-42 


49 


03 


,5 


19 


.87 





.01 


19 


,31 


0.01 


18 


,73 


0.01 


GC0444 


13 


26 


20. 


.63 


-42 


53 


46 





21 


.49 





.01 


20 


,80 


0.01 


20 


,15 


0.01 


GC0480 


13 


23 


50 


.49 


-43 


11 


43 


9 


20 


.96 





.03 


21 


,87 


0.03 


23 


.02 


0.02 


GC0481 


13 


23 


54. 


18 


-43 


21 


56 


6 


21 


.89 


0, 


.05 


22 


.81 


0.03 


24 


.10 


0.03 


GC0482 


13 


24 


03 


.98 


-43 


17 


20, 


3 


21 


73 





.03 


22 


,55 


0.02 


23 


.45 


0.03 


GC0483 


13 


24 


09 


.86 


-43 


22 


39. 





20. 


,89 





,05 


21 


,90 


0.03 


23 


.02 


0.03 


GC0484 


13 


24 


10 


.98 


-43 


12 


15 


,8 


20 


,53 





,02 


21 


.47 


0.01 


22 


58 


0.01 


GC0485 d 


13 


24 


13 


.76 


-43 


05 


33. 


,4 


19, 


,96 





,01 


20 


.83 


0.01 


21 


.84 


0.01 


GC0486 


13 


24 


19. 


.45 


-42 


54 


31 


,5 


18, 


,98 





,01 


19 


.83 


0.01 


20 


.77 


0.01 


GC0487 


13 


24 


20 


.60 


-43 


08 


07, 


9 


19 


.88 





.01 


20 


.83 


0.01 


22 


.01 


0.01 


GC0488 


13 


24 


27. 


.23 


-42 


57 


25, 


,7 


20 


,99 





03 


21 


,92 


0.03 


22 


.87 


0.02 


GC0489 


13 


24 


41 


.51 


-43 


12 


53 


2 


20, 


.30 





.03 


21 


,11 


0.02 


21 


,76 


0.01 


GC0490 


13 


24 


47. 


.40 


-43 


09 


59, 


,8 


20 


.20 





.01 


21 


,14 


0.01 


22 


,15 


0.01 


GC0491 


13 


24 


54. 


.98 


-43 


11 


36. 


9 


20 


.40 





.04 


21 


.41 


0.03 


22 


.46 


0.02 


GC0492 


13 


24 


56 


.61 


-43 


12 


23 


6 


20. 


.74 





.03 


21 


.42 


0.03 


22 


.00 


0.02 


GC0493 


13 


24 


57. 


.89 


-43 


07 


06 


3 


20. 


.11 





.03 


21 


.04 


0.02 


21 


,97 


0.01 


GC0494 


13 


24 


58. 


.30 


-43 


18 


33. 


6 


20. 


.20 





.02 


20 


.95 


0.02 


21 


,55 


0.01 


GC0495 


13 


24 


58. 


.96 


-43 


11 


28, 


,2 


19 


.98 





.02 


20 


.93 


0.01 


21 


.99 


0.01 


GC0496 


13 


24 


59 


.80 


-42 


57 


33 


,2 


19 


.68 





.01 


20 


.47 


0.02 


21 


.11 


0.01 


GC0497 


13 


25 


00 


.89 


-43 


21 


07 


,7 


20 


,37 





.01 


21 


.37 


0.01 


22 


,48 


0.01 


GC0498 


13 


25 


01 


,36 


-43 


05 


46 


,5 


20. 


,38 





.04 


21 


39 


0.01 


22 


.47 


0.01 


GC0499 


13 


25 


01 


.81 


-43 


03 


48, 


1 


20 


,43 


0, 


.02 


21 


.39 


0.02 


22 


.49 


0.03 


GC0500 


13 


25 


03 


.60 


-43 


01 


40, 


6 


20 


,42 


0, 


.03 


21 


.15 


0.01 


21 


,85 


0.03 


GC0501 


13 


25 


06 


.14 


-42 


46 


18, 


,8 


20 


,02 





.01 


20 


.91 


0.01 


21 


,97 


0.01 


GC0502 


13 


25 


06 


.68 


-43 


08 


58, 


,2 


20 


,98 





.06 


22 


.11 


0.05 


23 


,28 


0.04 


GC0503 


13 


25 


07 


36 


-43 


03 


23 


9 


19 


.08 


0, 


.01 


20 


,69 


0.01 


21 


18 


0.01 


GC0504 


13 


25 


08. 


.94 


-43 


08 


53 


,7 


20. 


,04 





.04 


20 


,70 


0.03 


21 


.29 


0.01 


GC0505 


13 


25 


12 


.69 


-43 


01 


55 


9 


19 


,80 





.04 


20 


,54 


0.03 


21 


.57 


0.03 


GC0506 


13 


25 


12 


.77 


-43 


10 


43 


,5 


20 


.42 


0, 


.04 


21 


,44 


0.04 


22 


.45 


0.02 


GC0507 


13 


25 


13 


.04 


-42 


55 


56, 


,2 


19 


,55 





02 


20. 


,32 


0.01 


20 


.94 


0.01 


GC0508 b 


13 


25 


13 


.70 


-42 


47 


36 


9 


21 


,27 





01 


22. 


,04 


0.01 


22 


.80 


0.01 


GC0509 


13 


25 


13 


.70 


-42 


40 


32 


4 


20. 


.73 





.01 


21 


,45 


0.01 


22 


,22 


0.01 


GC0510 


13 


25 


17 


68 


-42 


52 


56 


.5 


20 


.62 


0, 


.05 


21 


,44 


0.06 


22 


12 


0.02 


VjOOoll 


13 


25 


18 


.29 


-42 


47 


52 


.1 


20. 


.81 


0, 


.02 


21 


.66 


0.02 


22 


.62 


0.02 


LrOOoli 


13 


25 


18 


.92 


-42 


39 


03 


6 


19. 


.70 


0. 


.02 


20. 


.40 


0.02 


21 


.10 


0.01 


GO0513 


13 


25 


20 


.11 


-43 


22 


21, 


8 


20 


.33 


0, 


.01 


21 


.11 


0.01 


21 


.72 


0.01 


GC0514 


13 


25 


20 


.7-1 


-42 


41 


-17, 


1 


19 


,81 





01 


20 


,67 


0.01 


21 


.56 


0.01 


GC0515 


13 


25 


23 


.08 


-43 


06 


20 


9 


20. 


.65 





.02 


21 


,12 


0.02 


22 


11 


0.02 


GC0516 


13 


25 


23 


.31 


-42 


48 


17, 


.7 


19. 


.06 


0. 


.03 


19 


.88 


0.04 


20 


.93 


0.03 


GC0517 


13 


25 


24 


.75 


-43 


13 


37 


1 


20 


.40 





.01 


20 


.71 


0.04 


21 


.50 


0.02 


GC0518 


13 


25 


25 


.36 


-42 


49 


27, 


.5 


21 


.05 





.03 


21 


,75 


0.03 


22 


.15 


0.03 


GC0519 


13 


25 


26 


.09 


-43 


11 


22. 


.4 


18, 


,11 


0. 


.02 


18 


,81 


0.02 


19 


,50 


0.02 


GC0520 


13 


25 


27. 


18 


-42 


58 


37, 


9 


20 


,05 





03 


20 


,91 


0.02 


21 


.91 


0.01 


GC0521 


13 


25 


28. 


.81 


-42 


53 


15 


9 


19 


,82 





.01 


20 


,80 


0.02 


21 


.60 


0.01 


GC0522 


13 


25 


30 


.10 


-43 


22 


33. 


3 


20 


,44 





.02 


21 


.33 


0.01 


22 


.26 


0.01 


GC0523 b 


13 


25 


33 


.10 


-43 


09 


21 


6 


20 


,40 


0. 


.03 


21 


.16 


0.03 


21 


.76 


0.03 


GC0524 


13 


25 


33 


.94 


-42 


51 


39 


1 


20 


,85 





,02 


22 


,22 


0.01 


22 


.03 


0.02 


GC0525 


13 


25 


37 


.17 


-42 


43 


23 


.7 


20 


,77 





,04 


21 


.81 


0.04 


22 


,90 


0.02 


GC0526 b 


13 


25 


37. 


.63 


-42 


39 


32, 


2 


20 


,37 





,02 


21 


.11 


0.02 


21 


.93 


0.02 


GC0527 


13 


25 


39 


.41 


-42 


58 


24. 


1 


19 


,84 





,01 


20 


.44 


0.01 


21 


.08 


0.01 


GC0528 


13 


25 


48. 


.32 


-42 


55 


06 


9 


19, 


.91 





.02 


20 


.61 


0.01 


21 


,60 


0.01 




13 


9K 


so 


R.1 


-43 


OS 


09 


9 


90 


so 





09 


91 


m 


01 


91 


sn 


01 



597±52 


- 


543±10 


681±35 


- 


577±39 


487±102 


- 


363±32 


653±57 


- 


658±32 


807±65 


763±40 


526±87 


841±124 


719±31 


- 


666±140 


646±38 


805±60 


577±73 


499±51 


- 


247±49 


367±18 


410±38 


658±34 


657±32 


- 


705±120 


881±36 


- 


163±57 


- 


- 


401±55 


- 


466±96 


146±110 






503±75 






663±54 


593±40 


687±38 


396±84 






575±56 






505±51 


478±17 




706±63 






523±71 


440±15 


478±18 


490±113 






837±59 






700±48 






622±50 




649±42 


698±67 


597±11 


626±27 


615±81 




646±69 


225±35 






271±64 







646±61 

512±84 

599±105 

756±37 

498±22 

388±20 

566 ± 12 

579 ± 17 

559±35 

506±41 

568±22 

670±19 

713±28 

559±40 

676±42 

494 ± 17 

470 ± 33 

492 ± 25 

472±36 

600±35 

386±52 

564±71 

240±53 

670±41 

678±26 

578±20 

280±29 

460 ± 29 

326±71 

306±35 

670±45 

487±39 

664 ± 34 

844±54 

518±19 

591±29 

517 ± 12 

798±60 

589±63 

238 ± 25 

572±48 

160 ± 74 

573±49 

351±69 

527±61 

551±41 

507 ± 39 

262 ± 74 

489±19 

790+80 



503 ± 25 
553 ± 53 



619 ± 47 
635 ± 23 
437 ± 44 



478 ± 91 

441 ± 82 

581 ± 28 

250 ± 29 
179 ± 26 



519 ±45 
312 ± 22 



546± 9 

634± 26 

374± 30 

656± 27 

735±28 

726± 30 

690± 31 

524± 41 

354± 14 

648± 19 

858± 30 

163± 57 

417± 47 

146±110 

503± 75 

646± 24 

396± 84 

575± 56 

480± 16 

706± 63 

456± 10 

490±113 

837± 59 

700± 48 

637± 32 

603± 10 

632±52 a 

225± 35 

271± 64 

646±61 

512±84 

599±105 

756±37 

498±22 

388±20 

554±11 

576±16 

559±35 

506±41 

568±22 

670±19 

713±28 

559±40 

676±42 

508±15 

581±18 

478±21 

472±36 

600±35 

386±52 

564±71 

240±53 

670±41 

678±26 

578±20 

280±29 

461±27 

326±71 

306±35 

670±45 

487±39 

631±31 

844±54 

518±19 

591±29 

527±11 

798±60 

589±63 

243±19 

572±48 

176±24 

573±48 

351±69 

527±61 

551±41 

512±29 

301±19 

489±19 

790+80 
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TABLE 4 

Gaussian Fits for the Radial Velocity Distributions 



Group 


R 


mean 




sigma 


X 


red 




(arcmin) 


(km s- 1 ) 




(km s- 1 ) 




All GCs 


0-5 


516.1 ± 13 


.8 


i A c\ a I inn 

149.0 ± 1U.U 


1 


09 


An < 1 v 
All 


5-10 


ceo a n _l_ i n 
524.9 ± 10 


.3 


1 A O d _l_ T A 

14Z.0 ± / .4 





.87 


All GCs 


10-15 


525.3 ± 14. 


.3 


1 i I i n i 

152.1 ± 10.3 


1 


.15 


All GCs 


15-20 


461.0 ±24. 


.2 


177.7 ± 17.8 


1 


.80 


All GCs 


20-45 


520.4 ± 25. 


.2 


207.9 ± 18.6 


1 


16 


MP GCs 


0-5 


493.3 ±21. 


.7 


147.1 ± 16.1 





.88 


MP GCs 


5-10 


521.0 ±15 


.2 


141.4 ±11.0 


1 


.44 


MP GCs 


10-15 


506.3 ± 18. 


.8 


152.1 ± 13.5 


1 


.31 


MP GCs 


15-20 


446.0 ± 33. 


.5 


169.5 ± 24.4 


1 


.37 


MP GCs 


20-45 


475.1 ±31 


,5 


195.1 ± 23.6 


1 


.42 


MR GCs 


0-5 


529.6 ± 22 


.5 


166.8 ± 16.6 


1 


.27 


MR GCs 


5-10 


527.3 ± 14. 


.1 


144.0 ± 10.1 





.84 


MR GCs 


10-15 


550.0 ±21. 


.3 


148.1 ± 15.4 





.67 


MR GCs 


15-20 


474.4 ± 34. 


.4 


183.8 ±25.7 


1 


.33 


MR GCs 


20-45 


569.5 ± 38. 


.1 


204.1 ± 27.7 


1 


.74 



TABLE 5 

The Mass of NGC 5128 



GCs 




R 


N 


Vsys 


Ml 


e 


M p 


M r 


Mt 




(< 


iremin) 




(km s 


- 1 ) 


(km s~ 1 ) 


(deg. E of N) 


(xlO 11 M ) 


(xlO 11 M Q ) 


(xlO 11 M ) 


All GCs 


0- 


-5 


120 


514 ± 


14 


14 ±21 


175 ± 80 


1.900 ±0.663 


0.003 ± 0.009 


1.902 ±0.663 


All GCs 


0- 


■10 


317 


523 ± 


8 


33 ± 13 


194 ± 20 


4.520 ± 1.496 


0.028 ± 0.022 


4.548 ± 1.496 


All GCs 


0- 


•15 


432 


524 ± 


7 


36 ± 11 


190 ± 14 


6.257 ± 1.999 


0.050 ± 0.030 


6.307 ± 2.000 


All GCs 


0- 


■20 


487 


515 ± 


7 


36 ± 11 


189 ± 15 


9.260 ±2.881 


0.067 ± 0.041 


9.327 ±2.881 


All GCs 


0- 


■45 


549 


516 ± 


7 


33 ± 11 


185 ± 15 


17.45 ± 4.950 


0.119 ± 0.079 


17.56 ±4.949 


All GCs 


5- 


-10 


197 


528 ± 


10 


43 ± 16 


198 ± 19 


2.427 ±0.869 


0.047 ± 0.035 


2.472 ± 0.869 


All GCs 


5- 


■15 


312 


527 ± 


8 


44 ± 13 


192 ± 14 


4.128 ± 1.456 


0.075 ± 0.044 


4.203 ± 1.457 


All GCs 


5- 


■20 


367 


515 ± 


8 


44 ± 13 


191 ± 14 


5.832 ±2.038 


0.099 ± 0.059 


5.931 ± 2.038 


All GCs 


5- 


■45 


429 


517 ± 


8 


39 ± 12 


186 ± 14 


10.58 ±3.525 


0.166 ±0.102 


10.75 ± 3.527 


MP GCs 


5- 


-20 


180 


508 ± 


11 


46 ± 18 


206 ± 18 


4.879 ± 3.643 


0.109 ±0.085 


4.987 ± 3.644 


MP GCs 


5- 


■45 


216 


508 ± 


11 


31 ± 18 


192 ± 24 


9.942 ±7.153 


0.105 ±0.122 


10.05 ±7.154 


MR GCs 


5- 


■20 


184 


523 ± 


11 


50 ± 17 


174 ± 18 


6.157 ±2.564 


0.127 ±0.086 


6.284 ± 2.566 


MR GCs 


5- 


■45 


214 


525 ± 


11 


45 ± 16 


180 ± 18 


9.466 ± 3.766 


0.206 ± 0.146 


9.672 ±3.769 



TABLE 6 

The Kinematics of the Globular Cluster System of NGC 5128 



GCs 


R 


Raug 


N 


Vsys 


QR 


e 


@o ®major 




V.R/a Vp 




(arcmin) 


(arcmin) 




(km s~ 1 ) 


(km s- 1 ) 


(deg. E of N) 


(deg. E of N) 


(km s _1 ) 




All GCs 


0-5 


3.49 


120 


514 ±14 


14 ±21 


175 ± 80 


-40 


149 ±5 


0.09 ±0.14 


All GCs 


5-10 


7.31 


197 


528 ± 10 


43 ± 16 


198 ± 19 


-17 


147 ±3 


0.29 ±0.12 


All GCs 


10-15 


11.96 


115 


524 ± 14 


49 ±24 


184 ± 21 


-31 


151 ± 5 


0.32 ±0.17 


All GCs 


15-20 


17.32 


55 


446 ± 24 


43 ±38 


168 ± 43 


-47 


151 ± 5 


0.28 ±0.22 


All GCs 


20-45 


26.83 


62 


526 ± 22 


23 ±37 


55 ±87 


-160 


156 ±7 


0.15 ±0.25 


All GCs 


0-45 


10.82 


549 


517 ± 7 


33 ± 10 


185 ± 15 


-30 


150 ±2 


0.22 ±0.07 


MP GCs 


0-5 


3.78 


47 


496 ± 21 


43 ±29 


93 ±42 


-122 


144 ±8 


0.30 ±0.22 


MP GCs 


5-10 


7.41 


89 


524 ± 16 


46 ±24 


228 ± 27 


13 


147 ±6 


0.31 ±0.17 


MP GCs 


10-15 


11.97 


66 


511 ± 19 


47 ±33 


182 ± 30 


-33 


150 ±7 


0.31 ±0.23 


MP GCs 


15-20 


17.48 


25 


426 ± 31 


25 ±43 


137 ± 114 


-78 


128 ± 11 


0.19 ±0.42 


MP GCs 


20-45 


28.50 


36 


507 ± 28 


67 ±50 


45 ±37 


-170 


171 ± 14 


0.39 ±0.31 


MP GCs 


0-45 


11.90 


263 


506 ±9 


26 ± 15 


175 ± 28 


-40 


149 ±4 


0.17 ±0.09 


MR GCs 


0-5 


3.49 


57 


527 ± 22 


65 ±34 


244 ± 27 


29 


178 ± 11 


0.37 ±0.21 


MR GCs 


5-10 


7.23 


107 


527 ± 14 


51 ± 20 


178 ± 22 


-37 


149 ±5 


0.34 ±0.15 


MR GCs 


10-15 


11.96 


49 


542 ± 22 


46 ±36 


189 ± 34 


-26 


148 ±7 


0.31 ±0.26 


MR GCs 


15-45 


20.47 


54 


506 ± 24 


33 ±41 


168 ± 56 


-47 


162 ±9 


0.20 ±0.21 


MR GCs 


0-45 


10.20 


267 


526 ± 10 


43 ± 15 


196 ± 17 


-19 


156 ±4 


0.28 ±0.10 
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TABLE 7 

The Kinematics of Globular Clusters with Ages, Metallicities, and [a/FE] a 



Group N Ravg Vsys &o 

(arcmin) (km s _1 ) (km s -1 ) (deg. E of N) 



All GCs 


72 


5.7 


550±18 


25±26 


110±61 


GCs Age > 8 Gyr 


49 


5.8 


543±23 


27±36 


107±65 


GCs 5 < Age < 8 Gyr 


10 


4.9 


577±42 


53±78 


235±54 


GCs Age < 5 Gyr 


13 


5.7 


562±51 


58±92 


80±84 


GCs Age > 8 Gyr & [Z/H]> -1 


23 


5.9 


574±32 


17±50 


265±143 


GCs Age > 8 Gyr & [Z/H] < -1 


26 


5.7 


515±34 


55±48 


112±51 



a The ages, me t allicit ies and [a/Fe] used to eclassify these globular clusters are from 
IWoodlev etHI ([20101 1. 
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Catalog v p (km s _I ) 

Fig. 1. — We compare radial velocities measured from our thr ee surveys, LDSS-2 (t op panel), VIMOS (middle panel), and HYDRA 
(bottom panel) to the weighted radial velocity measurements in the lWoodlev et a l. (2007) catalog. The solid line is a 1:1 fit and the dashed 
line is the best least squares fit between the two sets of measurements. The least square fits are (slope, intercept)=(0.92, 60.2)££>g,g_2, 
(0.99, — 4.71)^/MOS) an d 0.81, 104. 1)hydra- We find good agreement in the LDSS-2 and VIMOS studies with the previously cataloged 
values. The HYDRA study has one significant outlier, which is GC0445 with cataloged velocity 504±45 km s _1 and our new measurement 
192 ± 39 km s —1 . GC0445 has been measured only one previous time in I Woodley et al. (2010). 
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Fig. 2. — A comparison of GCs that were measured more than once within each survey. The solid line is a 1:1 fit. LDSS-2 had 4 (solid 
circles), VIMOS had 1 (hexagon), and HYDRA had 17 multiple measurements (squares). The least squares fit for the Hydra data has a 
slope of 0.96 and an intercept of 13.1, indicating decent consistency of multiple measurements within the different HYDRA fields. 
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Fig. 3. — The radial velocity distribution of the entire GCS (open histogram), the metal-rich GC subpopulation (red dotted histogram), 
the metal-poor GC subpopulation (blue long dashed histogram), and the newly confirmed GCs from this study (hatched histogram). The 
entire GCS is fit with a Gaussian function using Rmix (mean = 516.7 ± 6.8 km s — 1 and sigma = 159.5 ± 4.8 km s — x ). See the electronic 
version for the color figure. 
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Fig. 4. — Radial velocity measurements for all metal-rich (red triangles) and metal-poor (blue squares) GCs as a function of projected 
galactocentric radius. See the electronic version for the color figure. 
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0-5' -rMP GCs 




200 400 600 800 200 400 600 800 200 400 600 800 

v p (km s" 1 ) 

Fig. 5. — Radial velocity histograms for the entire GCS (left panels), the metal-poor (middle panels) and metal-rich (right panels) GC 
subpopulations. Each subpopulation is binned radially from 0'-5' (top panels), 5'-10' (2nd panels), 10' -15' (3rd panels) ,15' -20' (4th panels), 
and 20'-45' (bottom panels). The best fits are listed in Table [4] 
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Fig. 6. — The C-Ti color of all GCs in NGC 5128 as a function of galactocentric radius. The dashed lines indicate the bounds of the GCs 
included in the kinematical and subsequent dynamical analysis. Excluded are the 4 reddest (GC0078, GC0408, GC0411, and GC0552) and 
the 3 bluest (GC0084, GC0282, GC0432) GCs. 
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Fig. 7. — Left: The luminosity function of all GCs in NGC 5128 with available Ti magnitude (N= 569) is fit by a Gaussian function with 
a mean of 19.44 ±0.04 mag and sigma = 1.03 ±0.03 mag. The luminosity functions for the metal-poor (blue) and metal-rich (red) GCs are 
also fit by Gaussian functions with means of 19.28 ± 0.06 mag and 19.36 ± 0.06 mag and sigmas of 0.90 ± 0.04 mag and 0.98 ± 0.04 mag, 
respectively. Right: The color-magnitude diagram of the entire GCS in NGC 5128 exhibits the color-bimodality of blue and red GCs. The 
dotted line represents the expected turnover magnitude. See the electronic version for the color figure. 
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Fig. 8. — The GCS of NGC 5128 is shown as their galactocentric radial position as a function of their azimuthal position measured in 
degrees East of North on the projected sky. The squares are the metal-poor GCs with a measured velocity, the triangles are metal-rich 
GCs with a measured velocity, and the solid circles are GCs that either have measured velocities but no color information, or GC s with no 
measured radial velocities. These latter objects were confirmed from resolved Hubble Space Telescope images (Harris ct al. 2006). Clearly, 
the outer regions (> 10') of NGC 5128 need to be searched for GCs to remove any potential spatial bias in this region. 
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Fig. 9. — The surface density profiles of the GCS is shown binned in circular annuli of equal numbers for the entire GCS (black circles), 
the MP GCs (blue squares), and the MR GCs (red triangles). These three populations were fit with a power law function between the 
regions of azimuthal coverage between 5'-20' indicated by the horizontal dashed lines. The best fits are labelled, and shown for the entire 
GCS (black line), MP GCs (long dashed blue line), and MR GCs (short dashed red line), respectively. The PNe power law fit between 
5'-20' is also shown (green dashed dotted line). See the electronic version for the color figure. 
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Fig. 10. — The total change in mass normalized by the number of GCs subsequently removed with decreasing v 2 R contribution in the 
Tracer Mass Estimator, is shown. We selectively remove 8 GCs with the largest v 2 R contribution indicated by the dashed line, which 
inflate the total mass estimate of NGC 5128. 
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Fig. 11. — The rotation axis and measured radial velocity are shown for the entire GC sample (top panel), the metal-poor sample (middle 
panel), and the metal-rich sample (bottom panel). The GCs have been fit with Equation [8| which is overplotted as the solid curve. The 
fitted parameters are listed in Table [6] 
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Fig. 12. — The v ays (top left panel), QR(top right panel), @o (bottom left panel), and cr Vp (bottom right panel) are plotted as a function 
of Ti magnitude. 



